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ABSTRACT 

We present a detailed analysis of the intrinsic UV absorption in the central HDFS target QSO 
J2233-606, based on a high-resolution, high S/N (~25 - 50) spectrum obtained with VLT/UVES. 
This spectrum samples the cluster of intrinsic absorption systems outflowing from the AGN at radial 
velocities v w —5000 - 3800 km s~^ in the key far-UV diagnostic lines - the lithium-like CNO doublets 
and Hi Lyman series. We fit the absorption troughs using a global model of all detected lines to 
solve for the independent velocity-dependent covering factors of the continuum and emission-line 
sources and ionic column densities. This reveals increasing covering factors in components with 
greater outflow velocity. Narrow substructure is revealed in the optical depth profiles, suggesting the 
relatively broad absorption is comprised of a series of multiple components. We perform velocity- 
dependent photoionization modeling, which allows a full solution to the C, N, and O abundances, as 
well as the velocity resolved ionization parameter and total column density. The absorbers are found 
to have supersolar abundances, with [C/H] and [0/H] «0.5 - 0.9, and [N/H] w 1.1 - 1.3, consistent 
with enhanced nitrogen production expected from secondary nucleosynthesis processes. Independent 
fits to each kinematic component give consistent results for the abundances. The lowest-ionization 
material in each of the strong absorbers is modeled with similar ionization parameters. Components 
of higher-ionization (indicated by stronger Ovi relative to Civ and Nv) are present at velocities 
just redward of each low-ionization absorber. We explore the implications of these results for the 
kinematic-geometric-ionization structure of the outflow. 
Subject headings: galaxies: quasars: absorption — galaxies: quasars: individual (QSO J2233-606) 



1. INTRODUCTION 

Active galactic nuclei (AGN) often show mass outflow 
as blueshifted absorption in their restframe UV and X- 
ray spectra. Many recent studies have explored the po- 
tential global effects of AGN outflows on all scales of their 
environment: from influe ncing the growth of the central 
supermassive black hole IjBlandford fc Begelman " 1999, 
12004' ) ■ the evolution of the host galaxy (Silk & Roes 
119981 IScannapJeco fc Oh 1 120041). the l uminosity func- 
tion of Quasars IjWvithe fc Loebl l2003() . and the mag- 
netization and chemical enrichment of galactic gas 
(jFurlanetto & Loeb "2001") and the intergalactic medium 
jCava lierc. Laoi. fc Mcnci 2002). It is unclear if the out- 
flows seen as UV/X-ray absorption are related to these 
global phenomena, or if they are connected to other AGN 
components showing evidence for bulk outflow motion 
such as the relativistic jets seen in radio observations 
and the emission-line gas. 

Measured ionic column densities provide the basis for 
physical interpretation of AGN outflows. Detailed UV 
spectral studies have shown measurements of these cru- 
cial parameters are often not straightforward. Analyses 
of absorption doublets and multiplets show the absorbers 
typically only partially occult the I j ackground emission 
sourc es (e.g., iHamann et al.l Il997t iBarlow fc Sargent! 
I1997J) . with potential complex scenarios involving dif- 
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ferent coverin g factors for d ifferen t background emis- 
sion sources l|Gangulv et all 119991: iGabel et all I2003D 
veloc ity-dependent covering factors (e.s.. lArav et alJ 
11999^. and possi bly inhomog eneous distributions of ab- 
sorbing material ijde Kool et al. 2002a; A rav et al. 20 05'). 
Detailed spectral studies of rich absorption line spectra 
at high spectral resolution and high signal-to-noise (S/N) 
are needed to deconvolve these effects from the observed 
absorption lines, and thereby determine the total column 
density, ionization, and chemical abundances of the out- 
flow via modeling. 

The quasar absorption-line spectrum of QSO J2233- 
606 (hereafter J2233) has been intensively observed 
and studied due to its location in the center of 
the Hubble Deep Field South (HDFS; Ferguson 
1998), its high redshif t (zp„,=2.238), an d relative 
brightness (B =17.5^ ISealev et a,].l 1199^ ISavagliol 
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1991 lOutram et all 119991: iProchaska fc Buries I 119991: 



Cristiani fc D'Odorico 1120001 IKim et al.ll200lD . In ad- 
dition to the rich intervening absorption spectrum, 
these observations have revealed seve r al "associated" 
absorption systems llSealev et al.l 119981: iSavaglio I Il998t 
lOutram et al]ll999t iPetitiean fc Srianandlll999D . These 
have been interpreted as outflowing material associated 
with the QSO in J2233 due to their partial line-of-sight 
coverage of the compact AGN emission, relatively broad 
widths, and proximi ty in redshift to the host g alaxy 
{zabs W2.198 - 2.21). IPetitiean fc Srianand I l)1999|) pre- 
sented an analysis of the (rest-frame) UV associated ab- 
sorption in J2233 using high-resolution optical spectra 
and supplementary long observations with the Hubble 
Space Telescope (HST) obtained as part of the HDFS 
program, which sampled lines spanning a large range of 
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ionization. Their photoionization modeling of one kine- 
matic component (zats =2.198) showed that in addition 
to a relatively low-ionization region giving the H I, C iv, 
and N v absorption, a high-ionization zone is needed to 
reproduce the strong Ne viii absorption. They found this 
component was consistent with having approximately so- 
lar metallicity, with qualitative evidence for the N/C 
abundance ratio enhanced relative to the solar value. 
They also found ion-dependent covering factors and sug- 
gested this may be due to a smaller covering factor of the 
emission-line r egion than the continuum source (see also 
lOiitram et a,l]ri9991. 

J2233 was observed for a total of 14 hours with 
VLT/UVES as part of a large program to s tudy inter- 
venin g systems in quasar absorption spectra (jKim et al.l 
I2001D . We present here a detailed analysis of the 
absorption-line properties and physical conditions in the 
associated absorption systems based on this spectrum. 
Our study is outlined as follows: in §2, we review the 
observational details and present the intrinsic absorp- 
tion spectrum from these observations. In §3, we make 
use of the broad spectral coverage which includes the 
CNO lithium-like doublets and H i Lyman series lines at 
high resolution and high S/N to solve for the individ- 
ual covering factors of the emission-line and continuum 
sources and the velocity-dependent ionic column densi- 
ties. In §4, we incorporate velocity-dependent photoion- 
ization modeling of each kinematic component to highly 
constrain the models and derive the full solution to the 
physical conditions (total column density and ionization) 
and chemical abundances of the outflow. In §5, we ex- 
plore the implications of this analysis for the outflow's 
geometric, kinematic, and ionization structure. 

2. VLT/UVES OBSERVATIONS AND THE INTRINSIC 
ABSORPTION IN J2233 

J2233 was observed with VLT/UVES for a total of 
50.4 ks between 1999 October 8-16 during the Com- 
missioning of UVES (we refer the reader to Cristiani & 
D'Odorico 2000 and Kim et al. 2001 for full details of the 
observations). These observations used the did 346x580 
(28.8 ks) and dic2 437x860 (21.6 ks) settings of UVES, 
giving full spectral coverage between 3060 - 10000 A. The 
data were processed with the ESO-maintain ed MIDAS 
ECHELLE/UVES package, as described in iKim et id] 
((2001). The reduced spectrum was resampled to 0.05 
A bins and normalized locally using a 5th order polyno- 
mial fit. The resulting final spectrum used in the analysis 
had resolving power ranging from R « 30,000 - 50,000 
and S/N between ~ 25 - 50 for the spectral regions of in- 
terest, with increasing values towards longer wavelength. 

The UVES spectrum covers the rest wavelengths of 
a host of key diagnostic far-UV lines for intrinsic ab- 
sorbers, including the CNO lithium-like doublets (O vi 
AA1032,1038, Nv AA1238,1242, Civ AA1548,1551) and 
Lyman series lines which are ubiquitous in intrinsic UV 
absorption spectra, as well as a number of less common 
lines that probe lower ionization (e.g., Cm A977, Siiv 
AA1393,1403, N in A991) and lower abundances (e.g., S vi 
AA933,944, Pv AA1118,1128). This is the first high- 
resolution spectrum of J2233 that spans all these lines 
with the same instrument and in the same epoch. The 
UVES spectrum is also of significantly higher quality 
than the previous observations: it has about twice the 



TABLE 1 

Intrinsic Absorption Components in QSO J2233-606 



Component 


Velocity'' 


FWHM'^ 


Lines'' 




km s^-*- 


km s^-"- 




1 


-4805 


56 


Ovi, Nv, Civ, Lya, Ly/3 


2 


-4137 


171 


Ovi, Nv, Civ, Lyo, Ly/3 


3 


-3971 


131 


Ovi, Nv, Civ, Lya'^ 


4 


-3830 


55 


Ovi 


5 


-3705 


50 


Ovi 


6 


-3634 


32 


Ovi 



'^ Centroid velocity and FWHM are for C IV A1548 for components 
1-3 and O VI A1032 for components 4 —6. 

^ Detected lines in each kinematic component in VLT/UVES 
spectrum. 

■^ Ly/3 contaminated with Lya forest. 



resolution and 4-5 times higher S/N at the C iv lines, 
twice the S/N at Nv and Lya, and six times better res- 
olution and 2 ~ 5 times higher S/N than the HST/STIS 
G430M grating observation of the O vi/Ly/3 spectral re- 
gion (see Table 1 in Petitjean & Srianand 1999). 

Figure 1 shows the associated absorption spectra for 
several key diagnostic lines in the UVES spectrum, plot- 
ted on velocity scale with respect to the restframe of the 
AGN. We also include the Neviii AA770,780 doublet ob- 
served with the Space Telescope Imaging Spectrograph 
(STIS) usi ng the E230M gr ating as part of the HDFS ob- 
servations (|Ferguson II1998() . Figure 1 shows three strong 
absorption components centered at Vr « —4800, —4150, 
and —3950 km s~^ (hereafter components 1, 2 and 3, 
respectively), identified as associated systems in earlier 
observations based on the partial line-of- sight covering 
factors inferred from the doublet lines llOutram et alJ 
119991 ISavagliollT99aHPe'titiean fc Srianand 1119991) . Thev 
are also distinguished from intervening absorbers in the 
high-resolution of the UVES spectrum by their relatively 
large widths, as seen for example by comparison with the 
Lya forest lines in Figure 1. These absorbers exhibit ab- 
sorption, in increasing strength, in the CiV, Nv, and 
Ovi doublets; they are not detectable in the lower ion- 
ization lines Cm or the Siiv doublet. Though the STIS 
spectrum has substantially lower S/N and there is poten- 
tial contamination from intervening absorption, NeViii 
also appears strong in components 1 and 3, and possibly 
component 2. Lya is present in all three components, 
and weaker than in the CNO doublets; Ly/3 is detected 
in component 1, with a potential weak notch present in 
component 2, while component 3 is contaminated with 
intervening Lya forest absorption. The higher order Ly- 
man lines are in a low S/N region of the UVES spec- 
trum, but there is no detectable absorption in these lines 
within the limits of the noise. There are three additional 
relatively weak systems seen in the Ovi doublet (and 
possibly Neviii) just redward of component 3 (compo- 
nents 4 - 6 at ^,. sa -3850 -3700, and -3625 km s'^ 
iPetitiean fc Sriana nd lll999l) . These systems are not de- 
tectable in absorption in any other lines. Table 1 summa- 
rizes the kinematic properties and lines present in each 
kinematic component. 
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Fig. 1. — Normalized absorption profiles from VLT/UVES spec- 
trum of QSO J2233-606. Each line is plotted as a function of ra- 
dial velocity with respect to the emission-Hne redshift of the quasar 
(z =2.238). The Nevill doublet from the HDFS E230M observa- 
tions is also included. Long wavelength doublet members and Ly/3 
are plotted with dashed lines, short wavelength doublet lines 
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Fig. 1. — and Ly a with solid lines. Regions plotted with dotted 
lines indicate contamination by intervening absorption or the other 
doublet member. Dotted vertical lines denote the three strong 
kinematic components analyzed in this study; tick marks show the 
three additional components seen only in O VI and Ne VIII. 



3. VELOCITY-DEPENDENT IONIC COLUMN DENSITIES 
AND COVERING FACTOR SOLUTIONS 

3.1. Fitting Formalism 

The observed intrinsic absorption troughs are a con- 
volution of hne-of-sight covering factor (C) and optical 
depth (r). These factors must be deconvolved in order 
to obtain ionic column densities for modeling the ioniza- 
tion and abundances of the outflows; additionally, the 
covering factors derived from this separation give im- 
portant constraints on the absorption-emission geome- 
try. With the two lines of an absorption doublet, single 
values for C an d r can be deriv ed at each vel ocity (e.g., 
iHamann et al.1 .1997: Barlow & Sargent Ill997f): however 
compl e x structure in th e absorber (e.g.. Ide Kool et aP 
l2002al lArav et S] I2QQ5D or the background emission 
.sources (e.g..lCa,ngulv et al .11799^ lOahel et a].ll2nn,'^ re- 
quires more than two lines to obtain the full solution. In 
particular, if continuum and line emission both undcrly 
the absorption features, the individual emission sources 
will in general have different covering factors, with im- 
plications for the doublet solution and relative values be- 



tween different lines (see §3 in lGabel et alJl2b05aj) . 

We applied a global fltting approach to determine the 
velocity-dependent covering factors and column densities 
of the outflow i n J223 3, following the method outlined 
in lGabel et alJ l)2005aj) . This approach allows for treat- 
ment of different covering factors for the continuum and 
emission-line sources by simultaneously fitting the ob- 
served absorption lines from multiple ions. The general 
absorption equation to solve for the normalized intensity 
of the j*'' fine is: 

/^=S4i?^(Gfe-^+l-Cf)], (1) 

where the i*'* individual emission source contributes a 
fraction i?^ = /^/Si[Jj^] to the total intrinsic intcn- 

sity and has cove ring factor Cj IjGangulv et al.l il999; 
iGabel et alJl2005al) . Here we have assumed that a single, 
uniform column density for each ion occults all emission 
sources. The effective covering factor for each line is the 
intensity-weighted combination of the individual cover- 
ing factors, 

c^ = s.cfi^^ (2) 

For components 1 - 3 in the J2233 outflow, there are 
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up to 8 detected lines available to constrain the absorp- 
tion parameters (fewer where lines are contaminated with 
other absorption): the C, N, O doublets, Lya, and Ly/3. 
We assumed all ions share the same covering factors of 
the continuum and emission line sources {Cc and Ci), 
and simultaneously fit the individual covering factors and 
Civ, Nv, Ovi, and Hi column densities. These param-o 
eters were solved in each velocity bin of the absorption^ 
profiles using equation 1 via the Levenberg-Marquardt m 
non-linear least-squares minimization technique (addi-g 
tional details are given in lGabcl ct al. 2005a) . We limited o 
our search of the covering factor parameter space to the 
physically meaningful range, 0< Ci <\. Regions of con- 
tamination with other absorption were omitted from the 
fitting. This included the far blue wing of component 1 
Nv A1242 and red wing of component 3 Nv A1238 due 
to blending of these lines, and obvious intervening Lya 
forest lines: at u w -4240 and -4140 km s~^ in Nv 
A1238 and v « -4030 km 8"^ in Ly/? (see Figure 1). 

A source of modeling uncertainty is the intrinsic 
flux levels of the continuum and BLR emission; the 
VLT/UVES spectrum is not flux-calibrated and thus the 
intrinsic shape of the AGN emission cannot be deter- 
mined from these data. Therefore, to determine the 
relative contributions. Re and Ri, we adopted results 
from previous observations of J2233. For O vi and Ly/3, 
we used the HDFS STIS G430M spectrum obtained in 
1998 October. For the longer wavelength hues, we used 
published results from earlier g round-based observat ions: 
Lya and Nv were from data in'Outr am et al.l (1999") ob- 
tained between 199 7 August - 1998 August and Civ 
from a spectrum in iSavaeliol (j[l998) obtained in 1997 
October. Based on these pubhshed results, our adopted 
contributions from the BLR for fitting were: Ri « 0.15 
for Ovi ; 0.3 for CiV ; 0.4 for Nv and Lya; and 0.1 
for Ly/3. A caveat is that since these observations were 
obtained at previous epochs (the UVES spectrum was 
obtained in 1999 October), there may have been vari- 
ability in the intrinsic emission that could affect our 
results. However, we believe this effect will be mini- 
mal since high-luminosity QSOs like J2233 do not typ- 
ically exhibit strong variability over timescales of a few 
years. For example, in a recent analysis of the w2500 
QSOs in the SDSS database having m ultiple observa- 
tions separated by more than 50 days, iWhilhite et al.l 
((^005) found only 12% exhibited significant variabil- 
ity. Additionally, HST/STIS observations of J2233 ob- 
tained a year apart had ide ntical fluxes (see Figure 2 in 
iPetitiean &: Srianandlll999j) . 

3.2. Results 

The results of the fits for the covering factors and op- 
tical depths are shown in Figures 2 and 3, respectively. 
The plotted error bars represent the formal la statis- 
tical errors in the best-fit parameters, corresponding to 
values giving Ax^ =1, and were computed from the di- 
agonal elements of the covariance matrix for the optimal 
fit (Bevington 1969); they do not include any uncertain- 
ties associated with the intrinsic flux levels. For cases 
where the covering factor solution is a boundary value 
(0,1), no covariance matrix elements are computed since 
it is not a minimum in the solution; for these cases, we 
estimated uncertainties by deriving solutions for models 
keeping the parameter fixed, and finding the value giving 
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Fig. 2. — Global-fit solution to the line-of-sight absorption cov- 
ering factors. Broad emission-line region covering factors are plot- 
ted with solid histograms, continuum covering factors with dotted 
histograms. Solutions were obtained by simultaneously fitting all 
detected lines in the UVES spectrum, assuming all ions have the 
same covering factors (see text). 



X^ =1 from the best-fit solution at the boundary value. 
For r >3, an absorption line can be considered to be 
saturated - for these large optical depths, changes in r 
affect the absorption depths less than the uncertainties 
in the flux levels. Thus, we adopt solutions giving r >3 
as lower limits of 3. This occurs over much of the core of 
the component 1 O vi profile, and over a narrower region 
of the core of Nv in this component. The fits in several 
bins in O VI in components 2 and 3 also reach these large 
values of r, mostly corresponding to the peaks in the r 
profiles in the other ions. In Figure 4, the absorption 
profile fits (red), derived by inserting the best-fit solu- 
tions into equation 1, are compared with the observed 
profiles. Most of the absorption in all lines is fit well by 
this model. In §4.5, we explore some of the minor dis- 
crepancies in the fit and their implications for our fitting 
assumptions and subsequent analysis. 

Figure 2 shows the continuum source covering factor 
is near unity for all velocity bins in components 1 and 
2. Component 3 only partially occults the continuum 
source, with Cc wO.75. The BLR emission is partially 
occulted in component 1, peaking at Ci «0.55 in the 
core and decreasing sharply in the wings of the profile. 
In components 2 and 3, the BLR coverage is low, with 
Ci =0 at better than the 2cr level at nearly all velocities. 
Results are summarized in Table 2. The overall covering 
factor solutions give compelling support for the assump- 
tions made in the calculations. Since each velocity bin 
provides an independent solution, the similarity in the in- 
dividual Cc solutions across each kinematic component 
supports the distinction of the continuum and BLR cov- 
erage as a major factor in the formation of the observed 
troughs. If this were instead an artificial distinction, it 
seems unlikely that the same result for Cc would be ob- 
tained over the numerous bins associated with each com- 
ponent, particularly since they exhibit a range of optical 
depths for each ion and optical depth ratios between ions. 
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TABLE 2 
Integrated Absorption Parameters 



1900 -4700 -4800 -4000 -380 

Velocity (km s"') 

Fig. 3. — Solutions to the intrinsic absorption optical depth 
profiles. Global-fit solutions are shown with black histograms. So- 
lutions to line pairs for each ion individually, assuming Cc from the 
global-fit solution, are shown as red histograms (see text). For the 
doublets, T for the red lines are shown. Tickmarks denote peaks in 
the narrow substructure seen in the solution to components 2 & 3. 



These results are also consistent with the general picture 
of the AGN emission geometry, in which the BLR is much 
more extended than the compact UV continuum source 
(further implications of these results are treated in §5). 
In Figure 3, interesting substructure in the optical 
depth profiles is revealed in components 2 and 3. The rel- 
atively broad features break down into a series of narrow 
subpeaks in r, with coincident kinematic structure seen 
in the different ions. We note some of the substructure in 
the optical depth profiles coincides with the fluctuations 
in the Ci solution in Figure 2. Thus, to test if the struc- 
ture in the r solutions is an artifact of the minimization 
fitting, we also fit the absorption constraining Ci —0 at 
all velocities. Resulting fits to the profiles are plotted in 
Figure 4 (dotted red histograms), showing similarly good 
fits over most of the absorption. The resulting r profiles 
of this fit show similar overall structure as the solutions 
shown in Figure 4 with Ci a free parameter, indicating 
this structure is likely real. Component 1 is seen to have 
velocity-dependent ionization structure; C iv and N v ex- 
hibit a narrow peak at w ~ —4810 km s^^, in contrast 
with O VI, which is strong (saturated) over a broader 



^a 


A^ovi'' 


A^NV 


Ac IV 


A^HI 


Afcin 


Cc'^ 


Ci'^ 


Id 


>13 


>7 


1.4 


1.1 


<0.1 


0.95 - 1 


0.55 


2 


19 


6.5 


1.4 


1.0 


<0.1 


0.9- 1 


0-0.2 


3 


15 


5.7 


1.2 


0.8 


<0.1 


0.65 - 0.75 


0-0.2 



^ Intrinsic absorption kinematic component (see Table 1). 

'' All column densities, Nion, are in units 10^"* cm~'^; O VI, N V, 
C IV, and H I were measured from the global-fit solutions (§3). 

•^ Derived continuum iCc) and emission-line (C;) covering factors at 
the core velocities of each component. 

'^ Heavy saturation in Ovi at i^ » —4800 4750; Nv saturated at 

V w -4820 - -4790. 



range of velocities extending redward to —4740 km s^^. 
H I lacks the sharp narrow peak in the core, exhibiting a 
relatively broad base of absorption over the profile. We 
integrated the optical depth profiles over velocity, and 
computed total ionic column densities associated with 
each kinematic component. Results are summarized in 
Table 2. 

Our measurements can be compared with those of 
Pctitiean & Srianand (1999), who present results only 
for component 1 . The H i measurements are similar 
(within Ri20%), however, our CiV column density is 
only about half their measured value, while Nv is at 
least two times greater, with saturation at some veloci- 
ties. Additionally, we find O vi is heavily saturated over 
much of the profile (see Figures 1 and 3) giving only 
a lower limit on its column density; this was not dis- 
cernible in the -ff^T/STIS G4 30M spectrum used in the 
iPetitiean fc Srianand I l)1999D study due to the hmited 
S/N. 

4. DERIVATION OF THE ABUNDANCES, IONIZATION, 
AND TOTAL COLUMN DENSITY IN THE OUTFLOW 

4.1. Modeling Methodology and Input Parameters 

To determine the physical conditions and abundances 
in the outflow in J2233, we compare the measured ionic 
column densities derived in §3 with predictions from 
models computed w ith the photoionization code Cloudy 
IjFerland et al.lll99^ . A model is specified by the spec- 
tral energy distribution (SED) of the ionizing continuum, 
the total hydrogen column density [Nh] and elemental 
abundances in the absorber, and the ionization param- 
eter [U = Q/AirR^nHc), which gives the ratio of the 
density of H-ionizing photons at the face of the absorber 
to the gas number density, uh- 

In components 1 - 3 in the J2233 outflow, model con- 
straints are available from measurements of C IV, N v, 
O VI, and Hi. To fully model these four ions, five pa- 
rameters must be specified for the adopted SED: U, Nh, 
and the abundances of C, N, and O relative to H, thus 
the solution for a single measurement of these ions is 
underconstrained. However, the high-resolution of the 
column density solutions derived in §3 can be utilized to 
give additional model constraints thereby allowing the 
full solution of the absorber properties. In the most 
straightforward model, the elemental abundances will be 
the same throughout the absorber; thus we model the col- 
umn densities measured in each velocity bin of an absorp- 
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Fig. 4. — Best-fit intrinsic absorption profiles for the eight fines 
solved with the global-fitting to t — C Model profiles (solid red 
histograms) were derived from the best-fit covering factors and 
column densities show in figures 2 and 3, and are plotted with the 
observed normalized profiles (black). 



tion component and link the solutions for all velocities by 
requiring the C, N, and O abundances to be the same at 
all velocities. Given the highly resolved absorption pro- 
files, this greatly increases the number of constraints for 
the models, which allows the fitting of the abundances 
and gives the velocity-dependent solution of the physical 
conditions {Nh and U) in individual absorption compo- 
nents. Another advantage of this method compared with 
modeling integrated column densities over a component 
is that regions of line saturation can be isolated. 

We adopted the SED used in Gabel ct al. (2005b) for 
our modeling, which consists of a broken power law 
(i^j/ ex i/~") in the ionizing continuum, with a = 1.4 
and 0.7 over energy intervals 0.0136 - 0.6 and 0.6 - 50 
keV, respectively. The ionizing continuum used in the 
calcul ations is a source of uncertainty; a study bv iLaorl 
(|1999) shows that the SED for quasars may be somewhat 
steeper than the one adopted for our modeling, which is 
typical for lower luminosity AGN. We have done exten- 
sive testing with different SEDs to explore its effect and 
find it does not significantly change our main results on 
the ionization and abundances in the outflow. 

We computed a grid of photoion ization models in U, 
Nh for solar elemental abundances IjGrevesse fc Anders I 



Fig. 4. — Model profiles for t — C solutions in which the BLR 
covering factor was fixed at C'l =0 are shown for components 2 and 
3 (dotted red), and are seen to be nearly identical as the solutions 
with Ci a free parameter. 



11989(1 . We then fit the absorbers by minimizing the x^ 
value in comparing the measured ionic column densities 
with the grid of model predictions. The U, Nh associ- 
ated with each of the n^ velocity bins are free parame- 
ters, with all velocity bins constrained to have the same 
CNO abundances relative to H. The abundances are in- 
troduced into the fit via a linear scale factor of the metal 
ionic column densities predicted by the solar metallicity 
grid of models. This is possible since all the absorbers 
considered here are optically thin to the EUV ionizing 
radiation (i.e., the Hen ionization edge), thus the total 
ionic column densities scale simply with the abundances 
of their parent element. In each velocity bin, there are 
measurements of Uion ionic column densities (up to four, 
depending on whether there is contamination or satura- 
tion present); the summation of these over velocity gives 
the total number of model constraints in each compo- 
nent. This gives 2 xrty+ 3 free parameters, to be mod- 
eled with Si=i,n„'T'ion;j Constraints. The best-fit solution 
is then the set of velocity-dependent U, Nh values and 
the C, N, and O abundances that minimize x^ from: 

2 ^ ^ y, r ^obs:j,i - Nmod-j{U,NH)fj .2 /gN 

where the measured column density of the j*^ ion in the 
j^th velocity bin Nobs-j,i is compared with the model col- 



20 

15 

10 

5 




VLT/UVES Observations of the AGN Outflow in HDFS Target QSO J2233-606 

_ . . , . . . . , . . . . ^^ 



Component 1 



- C 







5 


10 


15 




20 








c/c. 








20 


- 


IffU) }- 


'} 7 ^'" 


/^ 




.--■■■' — 


15 


- 


mM^ 


0^ 




:::n— 


- c - 


10 


r 


%^^^ 


:^- ■--'-' ''''' 






-_ 


5 












~ 















- 








5 




10 


15 


20 










c/c. 






20 


; 


- 


N 


//"/ 


^-^ 


_____o,oo-; 


15 


- 




/ 


/ / 


X^ 


- 


10 






// 


/ / /^ 




O- 1 


5 




/X 


/ / 




-^^^^ 


.-^ 


n 






== — 






- 



10 

N/Ns 



15 



20 



Fig. 5. — C, N, O abundance solutions from velocity-dependent 
photoionization modeling for component 1. The 2-D contour pro- 
jections of the 3-D reduced x^ volume are shown for each C-N-O 
plane. For each 2-D projection, the abundance of the third, perpen- 
dicular element was held fixed at its best-fit value. Solid contours 
show results for solutions to the column densities derived using the 
global-fitting to all lines; dotted contours show results for column 
densities derived individually for each ion pair, as described in the 
text. 



umn densities from the grid N„iod;j , fj is the scale factor 
relative to solar abundances, and aj^i the measurement 
uncertainty. For bins that are consistent with being sat- 
urated within noise limits (adopted to be r >3), we set 
the upper aji to a large value to effectively treat these 
as lower limits to the column densities. 

4.2. Abundance Solutions 

We have fit each of components 1-3 individually with 
this method, as well as the combined absorption from 
all components (i.e., constraining all kinematic compo- 
nents to have the same abundances). The analysis in- 
dicates super-solar abundances in all kinematic com- 
ponents, with generally consistent results for each ele- 
ment in each component. The solutions for each com- 
ponent individually are [C/H]=:0.6-0.7, [N/H] = l. 1-1.2, 
[O/H]=0.5-0.7, given as the log of the abundance ratio 
to hydrogen relative to solar values. Thus C and O are 
enhanced by a factor of « 3 - 5 relative to solar values, 
with additional overabundance of N consistent with Z^ 
scaling indicative o f enhanced secondary prod uction in 
massive stars (e.g., iHamann fc Ferland 1119991 and ref- 



FlG. 6. — Same as Figure 5, showing the results for components 
2+3 from modeling the global-fit column densities. 



erences there in). These results caii be co mpared with 
the study bv iPetitiean fc: Srianand I l)1999|) . who found 
qualitative evidence for enhanced [N/C], but no strong 
indication of supersolar C and O abundances based on 
their modeling. Figures 5 and 6 show the abundance 
solutions graphically for the independent fits to compo- 
nent 1 and components 2-1-3, respectively. This shows 
the 2-D (contour) projections of the 3-D x^ volume for 
the CNO abundances. The minimum reduced x^ value 
for the component 1 fit was 1.5, and for components 2-1-3 
«0.25. The low reduced x^ value for components 2 and 
3 may indicate the uncertainties derived in our column 
density fitting are overly conservative for these compo- 
nents. For each 2-D x'^ contour projection, the abun- 
dance of the other element was held fixed at its best-fit 
value. 

4.3. Velocity- Dependent U — Nh Solution 

Figure 7 shows the solutions to U and Nh as a func- 
tion of outflow velocity. The error bars represent values 
giving twice the minimum x^ value from equation 3 with 
the abundances fixed at their best-fit values. Compo- 
nent 1 is seen to have fairly constant ionization structure 
{log{U) ~ —1-2) over the blue-wing, and through the 
core of the component. The total hydrogen column den- 
sity increases from the outer wing, and levels off to a 
constant value over the central Aw «50 km s~^ in the 
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Fig. 7. — Velocity-dependent photoionization model solutions. 
The best-fit solutions are shown for the ionization parameter (mid- 
dle panel) and total H column density (bottom), with metal abun- 
dances shown in figures 5 and 6. Black, solid histograms show the 
solutions derived at the full resolution of UVES; red, dotted histro- 
grams show solutions where the data were binned by x2. Velocities 
of the narrow subpeaks identified in the optical depth profiles (Fig- 
ure 3) are indicated with tickmarks. The O VI doublet profiles are 
shown in the top panel for reference. 



core of the profile. In the red-wing, the total column and 
ionization increase sharply towards lower outflow veloc- 
ity. This can be seen in the absorption profiles, where 
O VI remains strong (saturated), while the other lines 
decrease. At w = —4770 km s~^, the ionization parame- 
ter has increased to log(C/)= —0.4 and Nh «10^^ cm~^, 
and by w = —4750 km s^^, there are only lower limits on 
these parameters due to the very weak absorption in the 
lower ionization lines. This may be due to the superpo- 
sition of two distinct physical components, or could rep- 
resent real physical velocity-dependent structure in the 
flow. For components 2 and 3, the physical conditions in 
the absorbers are not strongly velocity-dependent, with 

log{U) « —1.3 I.O. The total H column density does 

exhibit evidence for similar narrow substructure as the 
individual ionic column densities (§3), with Ri30% fluc- 
tuations. To test if this structure is an artiflcat of the 
solutions, in which strong fluctuations occur due to ef- 
fects of large optical depths (near the saturation limit), 
we have also computed models after binning the initial 
spectra by a factor of two. These results are plotted 
as grey, dotted histograms on Figure 7 for comparison, 
showing that much of the narrow structure is present in 
these smoothed solutions. 

We have generated synthetic absorption proflles for 
other lines observed in J2233 using the predicted ionic 
column densities from the photoionization model solu- 
tions and the covering factor solutions. These are shown 
in Figure 8 for key lines in the SST/STIS spectra and 
C III A977 in the VLT/UVES spectrum. Given the lim- 
ited S/N in the STIS E230M spectrum (Niv A765, Oiv 
A787, and Ne VIII A770) and the low-resolution of the O v 
A630 observation (G140L), it is difficult to derive strong 
conclusions from the HST observations; however, it is 
clear that the Ne viii is far underestimated by these mod- 
els, requiring a higher io nization component as demon - 
strated in the analysis bv lPetitiean fc Srianand 1 1)1999() . 
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Fig. 8. — Model absorption profiles for other observed lines in 
J2233. Synthetic profiles (heavy dashed lines) for lines not used 
in the photoionization modeling are compared with their observed 
profiles. They were derived using the predicted ionic column den- 
sities from the photoionization model solutions and the covering 
factor solution. Lines from the HST/STIS E230M spectrum (N IV, 
O IV, and Ne vill) were smoothed over 5 bins; the model O v pro- 
file was binned to approximate the G140L resolution. The Cm 
spectrum is from the VLT/UVES observation. The VLT/UVES 
spectrum of O VI is shown at the top for reference. 



A similar situation was found in the case of QSO PG 
0946+301 from a detailed analysi s of the highly con- 
strained BAL bv lArav et alJ lj2001a|l - an additional high- 
ionization absorber coincident in velocity with the other 
UV lines was needed to reproduce the observed Neviii 
absorption. The other lines in J2233 are consistent with 
the observations, with the possible exception of N IV A765 
in component 2 which may be overestimated. The model 
predictions for Cm A977 are also seen to be consistent 
with the observations. Additionally, absorption from the 
Si IV A1393,1403 doublet is predicted to be negligible in 
these models and consistent with the tight upper limits 
on the Si IV column densities from the UVES spectrum 
(see Figure 1). We note that for the parameters derived 
here, these UV absorption systems will produce negligi- 
ble X-ray warm absorber features. 

4.4. High- Ionization Components 4^6 
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Here we derive constraints on components 4-6. With 
only Ovi detectable in the UVES spectrum, there are 
not enough constraints to solve for the individual cover- 
ing factors. However, comparison of the depths of the 
doublet members reveal the weak component 4 just red- 
ward of component 3 has low covering factor. Though 
uncertainties from the doublet s olution are rather la rge 
for weak features like this fsee iGabel et alJl2005aj) . it 
requires that the continuum source is only partially cov- 
ered (i.e., even if the BLR is fully unocculted), with Cc 
at least as low as the value derived for component 3. 
The effective covering factors for components 5 and 6 are 
larger, consistent with C ^ 0.8 - 1. We note that these 
solutions should be taken with caution due to possible 
contamination with the Lya forest, and with no other 
lines for confirmation. 

To constrain the physical conditions in these absorbers, 
we mea sured O VI column de nsities using the doublet so- 
lution IjHamann et al J 119971) . and derived upper limits 
on Hi, Nv, and Civ. This gives lower limits on U, Nh 
via photoionization modeling. The ionization in all three 
components is substantially higher than components 1 - 
3. Adopting the abundances derived above for compo- 
nents 1-3 implies \og{U) > —0.25, 0, and —0.4 and 
log(iV//) > 17.7, 18.4, 17 for components 4-6 respec- 
tively, but the values could be substantially greater. Up- 
per limits on U, Nh come from the H i upper limit, giving 
\og{U) < 0.5, log{NH) < 19.7 for each component. As- 
suming the features in the Ne viii STIS spectrum are real 
for these components, the implied column densities are 
consistent with the upper range. 

4.5. Implications of More Complex Absorption 
Scenarios 

As described in §3.1, our solution to the absorption pa- 
rameters {t,C) involved some simplifying assumptions. 
Specifically, we assumed: (1) the covering factors are 
identical for all ions (separately for the continuum and 
BLR sources), and (2) the absorption column density for 
each ionic species is uniformly distributed over the oc- 
culted fraction of the background emission (i.e., it is a 
perfect step-function in our sightline to the emission re- 
gions). We explore the implications of these assumptions 
in more detail below. 

Ion-Dependent Covering Factors 

The assumption that all ions have the same covering 
factors was invoked in our global fitting to increase the 
number of fitting constraints so that the individual con- 
tinuum and BLR covering factors could be solved. The 
results from this fitting matched the observed profiles 
well overall (see Figure 4), indicating the assumption 
is not too far off. However, there are some discrepant 
regions in the fit which may signify a more complex 
absorption-emission geometry in the J2233 outflow; pri- 
marily, the absorption strength in the C iv blue (red) 
doublet member is systematically underestimated (over- 
estimated) over the core region of component 1 and in 
some regions of component 2. 

To explore this further, we have computed column den- 
sities and covering factors for each doublet/Lya-Ly/3 pair 
individually using insight from our global fitting. Specifi- 
cally, we fixed the continuum covering factor at the aver- 
age value for each kinematic component from the global 



fitting {Cc = 1, 0.92, 0.75 for components 1, 2, and 3 
respectively; see Figure 2), and then solved for r and Ci 
separately for each ion line pair by applying x^ analysis 
to equation 1. We believe this is the best approach to 
explore ion-dependent covering factors since the BLR is 
much larger than the UV continuum source, thus differ- 
ences in covering factors between ions would most likely 
be due to affects of BLR coverage. For la uncertainties 
in these parameters, we adopted the maximum offsets 
from the best-fit solutions which give x^ =2, since there 
are two lines being fit. Solutions can only be obtained 
at velocities where both members of the line pair are 
uncontaminated with other absorption. Bins that are 
saturated (we adopted r > 3) have only lower limits on 
their column densities. In Figure 3, we plot the column 
density profiles derived from this analysis. The primary 
differences of these solutions are that the C iv column 
density is somewhat larger in component 1 (and equiv- 
alently Ci is lower), by up to «50%, and the associated 
1(7 uncertainties for all ions are generally larger than in 
the global-fit solutions. Figure 3 shows that most of the 
narrow substructure in the optical depth profiles in com- 
ponents 2 and 3 described in §3.2 is still present in the 
individual ion solutions. 

To test the potential affect of ion-dependent covering 
factors on our photoionization modeling results, we did 
the analysis described in §4.1 using the column densities 
derived individually for each ion in component 1 (most 
of components 2 and 3 have no solutions for either N V or 
H I due to contamination; see Figure 3). The results for 
the abundances are presented in Figure 5, showing the 
X^ values as dotted contours. The carbon and nitrogen 
abundances derived from this model are about a factor 
of 1.5 times greater than the solutions using the global- 
fit column densities, whereas the oxygen abundance is 
similar. We note that since O vi is saturated over a large 
portion of the component 1 profile and thus has only 
lower limits on the column density, this may lead to an 
underestimate of the oxygen abundance. 

Inhomogeneous Absorption 

Another assumption we have made in our t — C 
measurements is that the absorber is a perfectly ho- 
mogeneous slab, with uniform column densities span- 
ning the occulted portion of the background emission. 
More complicated distributions of outflow material, with 
a continous range of optical depths distributed with 
different covering factors (i.e., an inhomogeneous ab- 
sorber), were first explored by de Kool et al. (2001). 
More recently, [Ara v et al. (2005, hereafter AK05) and 
iSabra fc HamanrTI (2006, hereafter SH06) presented de- 
tailed analyses of inhomogeneous absorber (lA) distri- 
butions and their implications for observed absorption 
troughs. Results from these studies indicate inhomoge- 
neous absorption will not have a large affect on our mod- 
eling results and conclusions for the outflow in J2233, as 
described below. 

In AK05, it was shown that for lines from the same 
ionic level that have different optical depths (e.g., dou- 
blet line pairs or Lyman series lines), cases where the 
observed residual fluxes in the two line troughs are finite 
and equal in high S/N spectra (i.e., non-black satura- 
tion in the traditional partial coverage models), cannot 
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be explained straightforwardly by lA models. Instead, 
they require a sharp edge in the distribution of absorber 
material, consistent with the homogeneous partial cover- 
age model. Specifically, using simple functional forms for 
the absorber material (power-law and Gaussian distribu- 
tions), AK05 showed that "shallow" distributions (i.e., 
power-laws with small spectral indices or Gaussians with 
relatively large b parameters) cannot produce finite but 
equal residual fluxes for any conditions, and steep distri- 
butions require unrealistically high values for peak opti- 
cal depths (see their Figures 2 and 3). The more general 
calculations of SH06, which allowed non-zero minimum 
bounds on the optical depth distribution, alleviated the 
peak r requirements somewhat for cases of low residual 
fluxes. However, even for these calculations, there is no 
way to obtain an absorption profile for a single line pair 
that gives equal residual fluxes at some velocities but 
different fluxes at others unless the functional form of 
the distribution of material is velocity-dependent across 
the component (e.g., see the models with C/ =0.9 for 
the HPC models in Table 1 of SH06). This is precisely 
the case observed in the troughs of J2233, for example 
in O VI and N V in component 1. Such a situation 
seems physically implausible. Thus we conclude that the 
lA models as presented in the detailed studies of AK05 
and SH06 do not provide a natural interpretation for the 
J2233 outflow. 

Finally, we note SH06 have shown that even for cases 
where inhomogeneous absorption does apply, the mean 
optical depths derived by assuming homogeneous par- 
tial coverage are very similar to the mean values for 
lA distributions for most cases (they diverge somewhat 
for strongly peaked distributions). Thus, the solutions 
from the homogeneous partial covering analysis are gen- 
erally representative of the inhomogenous distributions 
(for cases that are not saturated), once the mean opti- 
cal depths are accounted for by averaging over the unoc- 
culted and occulted regions. It has yet to be shown, how- 
ever, the implications for photoionization models which 
are done assuming a slab of uniform column density. 

5. IMPLICATIONS FOR THE GLOBAL NATURE OF THE 
J2233 AGN OUTFLOW 

Here we synthesize the above results on the kinemat- 
ics, geometry, ionization, and abundances to explore the 
global nature of the outflow in J2233. 

5.1. Chemical Enrichment 

Determining abundances in quasar outflow systems 
serves as a probe of the chemical evolution of the uni- 
verse, with implications for the timing and mechanisms 
for nucleosynthesis and distribution of metals on galac- 
tic and intergalactic scales. Analysis of intrinsic narrow 
absorption-line systems provides the best means for prob- 
ing the abundances of QSO environments: BAL lines are 
typically heavily saturated and blended preventing de- 
termination of ionic column densities, whereas observed 
emission-lines sample the integrated flux over a large spa- 
tial region, and thus likely sample a large range of phys- 
ical conditions. Evidence for near-solar or super-solar 
abundances has been found in associated absorbers in 
se veral other hi gh rcdshift (z ^ 2-4) QSOs (see references 
in lHamann fc Ferland 1999^ iD'Odorico et al.i 2004) . al- 
though this is often not definitive due to limited quality 



of spectra. Additionally, linking these results to QSO 
outfiows is uncertain in many cases since they lack other 
properties that would indicate the absorbers are intrinsic 
to the AGN. 

As shown in §4.2 and Figures 5 and 6, all three outflow 
components in J2233 with available fitting constraints 
have supersolar abundances, with enhanced N enrich- 
ment: Z « 3-5 Zq for C and O, and N/H scaling roughly 
as Z^. The similarity in results derived independently in 
each kinematic component gives strong support for this 
finding. Additionally, the enhanced N abundance rela- 
tive to C and O supports the general finding of high abun- 
dances in the absorbers, since secondary CNO nucleosyn- 
thesis becomes important after the system has reached 
near solar metallicity, based on ob servations a nd chem- 
ical enrichment simulat ions (e.g., iHamann fc Ferlandl 
Il993l iVUa-Costas "fcEd munds 199% The advanced 
chemical processing implied by our detailed fitting in- 
dicates the host galaxy of J2233 underwent vigorous 
star formation at an early epoch, before z =2.2. Given 
the finding of high metallicity in associated absorbers 
in other bright, distant QSOs, this supports the evolu- 
tionary scenario in which massive galactic nuclei undergo 
rapid, extensive star formation at early times, with most 
chemical enrichment occurring before at least z > 2-3 
(Hamann & Ferland 1999, and references therein). A 
subsequent or coeval QSO phase may then distribute the 
enriched gas to the host galaxy, and possibly IGM, via 
an outfiow. The finding of distances of AGN outflows 
from the nucleus on kpc scales or greater l)de Kool et alJ 
2001; Hamann et al. 2001) supports this general model 
that QSO winds can seed their environments with met- 
als. 

5.2. Kinematic- Geometric- Ionization Structure 

The detailed results of our absorption fitting and mod- 
eling reveal correlations in the kinematic, geometric, and 
ionization structure in the J2233 outflow. For example, 
our analysis in §3 (Figure 3) shows the covering factors of 
the individual kinematic components with low-ionization 
absorption are correlated with outflow velocity. To sum- 
marize, the highest velocity absorber, component 1, fully 
occults the continuum source and covers w 60% of the 
BLR, component 2 is consistent with full continuum cov- 
erage but little or no BLR coverage, while component 3 
only partially occults the continuum source («70%). 

One interpretation of this trend is that the absorp- 
tion region size increases with outflow velocity, at least 
as projected on the sky against the background AGN 
emission (i.e., in the transverse dimension, Xt)- //the 
absorbers in our sightline to the nuclear emission have 
a uniform density filling their entire volumes (i.e., filling 
factor / =1), e.g., if each component is a single "cloud" 
or continuous flow region, then the expression relating 
the H column density to the volume density and radial 
thickness of the absorber is Nh — n„/S.R. Since the 
integrated Nh for the low-ionization regions are simi- 
lar for components 2 and 3 and somewhat greater for 
component 1, if the radial dimensions of the absorbers 
follow the same trend as the transverse sizes in the sce- 
nario described above (e.g., for a cloud-like geometry), 
then the volume densities decrease with outflow velocity. 
Furthermore, the similarity in ionization parameter for 
the low-ionization material in each component (§4, Fig- 
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ure 7) implies their uhR^ products are similar via the 
expression for U given in §4.1. Combining these results 
would imply a correlation between outflow velocity and 
R for the components in J2233 if all of the above as- 
sumptions apply. Qualitatively, this result is consistent 
with the predictions from dynamical models which have 
continuous acceleration of a radial outflow. Another in- 
terpretation of the C — Vout correlation is that outflow 
regions with higher velocity have moved further across 
our sightline to the emission sources. 

5.2.1. Geometric Constraints 

The distance of the absorber from the central ionizing 
source {R) has important implications for assessing the 
total energy and mass flux associated with the outflow. 
The most direct and accurate way of determining R is 
by measurement of the volume density in the absorbers, 
n„, from excited level ab sorption (e.g., Ide Kool et al.l 
l20f)ll I2002bic]: ICxabel et alJl2005bt) . but the presence of 
these diagnostic lines in AGN outflows is rare. They are 
not detectable in the spectrum of J2233 within the noise 
limits; this is due to the relatively low- ionization of the 
species giving excited level UV absorption (e.g., Sin* 
A1260, Cii* A1336, Cm* A1175) compared with the ion- 
ization state of the outflow. However limits on R and the 
geometrical structure of the outflow can still be explored 
based on indirect arguments involving other measurable 
quantities. 

For example, if the filling factors of the absorbers are 
unity as described above, then combining the expres- 
sio ns for Nh and the ionization parameter (§4.1) follow- 
ing iTurner et alJ (p993) gives the geometrical parame- 
ters R and Ai? in terms of the parameters derived from 
photoionization models, the H-ionizing photon flux, and 
constants: 

Ai? NhUAt:c 



i?2 



Q 



The photon flux can be estimated from the luminosity 
distance to J2233 (20 Gpc for Ho = 71 km s."^ Mpc"!, 
rim =0.27, r^A =0.73; Spcrgcl et al. 2003) and integra- 
tion of our model SED over the ionizing continuum, giv- 
ing Q =3x10" s'^. 

We can use the covering factors derived for the BLR 
(§3) and the size of the BLR to constrain the projected 
transverse size {Xt) of the absorption regions in our 
sightline to the background AGN emission. Usin g the 
relation {Rblr oc L^^,; e.g.. iWandel et aUlTOfl^. this 
gives Xt > 1 pc for component 1 [Ci «0.5), and smaller 
for components 2 and 3 (C; «0). If the radial sizes of 
the absorption regions are similar to these limits on the 
transverse sizes, Ai? « Xt, consistent with a cloud-like 
geometry, then equation 4 implies i? « 50 kpc if the flU- 
ing factor of the absorbing medium is unity. This large 
distance seems improbable, requiring the presence of ex- 
tremely small structure relative to its distance from the 
nucleus, with just the right line-of-sight geometry to par- 
tially cover the AGN emission. However, similar results 
have been implied in analysis of other QSO outflows (e.g., 
iHamann et alJl2001|) . Alternatively, if the outflow has a 
thin shell geometry (Ai? ^ Xt), it could be consider- 
ab ly closer in to t he nucleus; this model was adopted 
bv lSteenbrugge et al. (.2005,), for example, to explain the 
X-ray absorption in NGC 5548. Another alternative is 



that the volume filling factor of the absorption material 
in our line-of-sight is small (i.e., it is highly clumpy). In 
this case, the effective radial thickness of the absorber is 
proportional to f^/^l\R, placing the outfiow closer to the 
nucleus for a given Ai? : Xt ratio. We note that very 
low volume filling factors appear to be a generic feature 
of QSO outflows seen in the UV, and its implications and 
requirements on dynamical models was recognized early 
on in the study of the se systems fe.g.. IWcvmann et all 
IT985llAra7erai]IT99l . 

5.2.2. Narrow Kinematic Substructure and 
Velocity- Dependent Ionization 

An interesting feature of the outfiow in J2233 is the 
prevalence of narrow kinematic structure seen in the ab- 
sorbers. This is directly apparent in the clustering of dis- 
tinct kinematic components 2-6, which combine at ad- 
jacent velocities to give continuous absorption over Aw > 
500 km s~^. Also, the detailed fitting and modeling in 
§3 and 4 revealed the relatively broad component 2 and 
3 absorbers are comprised of a series of narrow peaks in 
column density (Aw w 25 - 50 km s^^). Further, compo- 
nent 1 shows strong velocity-dependent ionization differ- 
ences: the red-wing is of significantly higher ionization, 
with saturated O vi absorption and weak or negligible 
absorption from lower-ionization species. This may be 
interpreted as the superposition of two physically distinct 
components, partially overlapping in velocity space. 

Similar clustering in velocity of multiple, relatively 
narrow components is commonly seen in QSO outflows. 
For example, many BALQSOs show narrow substructure 
in some of the weaker lines coinciding in velocity with 
or adjacent to the broad tro ughs ('e. g..lVoi t et alJll 993t 
lArav et al.lll999Hde Kool et al. 2001; Arav et al.l2001b(l . 
In a recent m ulti-epoch study of QSO HS 1603+3820, 
iMisawa et alJ 1)2005') found a cluster of partially overlap- 
ping narrow components had become deeper and more 
blended, appearing almost as a BAL trough. These re- 
sults indicate we are seeing a number of distinct absorp- 
tion regions in our line-of-sight to the AGN emission, 
either different stream lines associated with continuous 
outflows or distinct clouds. Perhaps there is a continuous 
distribution of absorption widths associated with AGN 
outflows, ranging from the narrow components seen in 
J2233 to the classical BAL troughs, with the observed 
width related to the details of the outflow mechanism or 
orientation of our sightline to the nucleus. The cluster- 
ing in velocity suggests the distinct regions are related, 
possibly determined by the details of their origin or driv- 
ing mechanism. It is interesting that the low-ionization 
component 1 and 3 absorbers in J2233 both have higher 
ionization absorption at adjacent, lower outflow veloci- 
ties. Similar structure was seen in the QSO 2359-1241 
absorber, and was taken as evidence fo r ionizat ion strat- 
iflcation in a disk wind bv ' Arav et alJ l|2001bD . In con- 
trast, Voit et al. ( 1993) found some objects show the op- 
posite trend, with low-ionization absorption coinciding 
in velocity with the low-velocity region of BAL troughs. 

6. SUMMARY 

We have presented a detailed analysis of the intrinsic 
UV absorption in the central HDFS target, QSO J2233- 
606. Our analysis is based on a high-resolution {R k, 
30,000 - 50,000), high S/N (- 25 - 50) spectrum obtained 
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with VLT/UVES. This spectrum samples the cluster of 
intrinsic absorption systems outflowing from the AGN at 
radial velocities v « —5000 - 3800 km s~^ in the lithium- 
like CNO doublets and Hi Lyman series lines, which 
serve as key diagnostics to determine the covering fac- 
tors, column densities, ionization state, and abundances 
of the outflow. 

We fit the absorption troughs as a function of radial 
velocity using a geometric model which assumes inde- 
pendent line-of-sight covering factors for the individual 
background emission sources (continuum and emission- 
line) and the same covering factors for all ions. The 
solution reveals: 

1. There is a trend of increasing covering factor for 
components with greater outflow velocity: compo- 
nent 1 at u w —4950 km s'-^ occults the entire 
continuum source and Ri55% of the BLR; compo- 
nent 2 at w w —4150 km s^^ covers the continuum 
source, but little or none of the BLR; component 
3 at w ss —3950 km s^^ only partially covers the 
continuum source («70%). 

2. Narrow kinematic substructure in the ionic column 
densities is prevalent. The red-wing of component 
1 is highly ionized, having saturated O vi and weak 
absorption from lower ionization species compared 
with the core of this component. Components 2 
and 3 break down into a series of subpeaks in r, 
seen consistently in the solutions to different ions. 

We performed detailed photoionization modeling to de- 
termine the metal abundances and physical conditions 
in the absorbers. Ionic column densities were simulta- 
neously fit for all velocities of each component by con- 



straining the abundances in each velocity bin to be the 
same. This technique greatly increased the number of 
modeling constraints, allowing a fit to the individual ele- 
mental abundances of C, N, and O relative to H, as well 
as the velocity-dependent ionization parameter and total 
column density. Our analysis shows: 

1. All absorption components have supersolar metal 
abundances, with [C/H] and [0/H] «G.5 - 0.9, and 
[N/H] « 1.2. This is consistent with the Z^ scaling 
of N predicted from secondary CNO nucleosynthe- 
sis. Consistent results are found for independent 
fits to each component individually. 

2. The lowest ionization material in each of compo- 
nents 1-3 has similar ionization (log(C/) « —1.2). 
Higher ionization gas is required at the same veloc- 
ities to account for the Ne Viii absorp tion, as shown 
also bv lPetitiean fc Srianand 1 1)1999() . The redwing 
of component 1, and components 4-6 which lie 
just redward of of component 3, have ionization 
parameters at least 5-10 times greater than the 
low-ionization component 1-3 gas. 
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